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We consider quantum dots with a parabolic confining potential. The qualitative features of 
such mesoscopic systems as functions of the total number of electrons N and their total angular 
momentum J, e.g. magic numbers, overall symmetries etc., are derived solely from combinatoric 
principles. The key is one simple hypothesis about such quantum dots yielding a basis of states 
(different from the usual single electron states one starts with) which is extremely easy to handle. 
Within this basis all qualitative features are already present without the need of any perturbation 
theory. 

PACS numbers: 73.20.Dx, 73.23.-b, 73.61.-r, 71.10.-w 
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I. INTRODUCTION 

Quantum dots with a parabolically confined potential 
are a particularly simple species of a phenomenon en- 
joying increasing interest The basic novelty with 
such systems is their being mesoscopic, i.e. the fact that 
they live on the edge between classical and quantum 
physics. Quantum dots are sometimes called artificial 
atoms, which illuminates their mesoscopic nature. In- 
deed, although much larger than a real atom, the angular 
momentum j of the electrons is quantized, and hence, via 
the classical relation j oc r 2 , its orbital radius. 

It has been known for some time [|]||] that the to- 
tal potential energy of a (parabolic confining) quantum 
dot with N electrons exhibits local minima for certain 
amounts of the total angular momentum J, called magic 
numbers. In fact, the local minima occur when 
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The quantum nature of the system becomes even more 
apparent when we compare the behaviour for N < 5 elec- 
trons with the one for N > 6 electrons. Classically, we 
expect that the lowest energy configuration for N < 5 
electrons is a regular iV-gon, while for N > 6 electrons 
it should be a regular (N — l)-gon with one electron 
in the center. However, whenever k = (N — l)k' with 
k' G Z + , the configuration at such special magic numbers 
is much less symmetric than the classical counterpart. In 
fact, it behaves much more like a Laughlin-type quantum 
droplet. 

We will demonstrate in this letter that the appearance 
of magic numbers as well as the loss of symmetry for the 
special magic numbers satisfying 



J(N, k'(N - 1)) = N(N - 1) 



k' e Z+ , (2) 



can be explained by pure combinatorics. To do so, we 
need a bit of preparation. 

We assume that the external magnetic field is strong 
enough to completely polarize all electrons. Then, only 
one electron can occupy an eigenstate of the angular mo- 
mentum operator with eigenvalue j € Z + . The system 
characterized by the total number of electrons N and the 
total angular momentum J may now reside in any quan- 
tum mechanical state of the form \j\ , . . . jjf] such that 

< ji < h < ■ ■ ■ < ]n and YliLi Ji = J > where the j. t 
denote the angular orbitals in which the electrons may 
stay. Of course, the N electrons are identical particles 
meaning that each angular momentum eigenstate is a 
one-particle state antisymmetrized over all the electrons. 
To count the number of possible configurations one sim- 
ply has to expand the combinatoric partition function 



Z(y,q) = l[(l + y q n ) 
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The product includes the term n = since one electron 
may reside in the j ' — state. The reader might convince 
themselves that cn,j can be quite large even for very 
small electron numbers and medium sized total angular 
momenta. The true state is a complicated superposi- 
tion of these configurations which depends on the precise 
interaction of the electrons giving each configuration a 
different weight. The problem is that each configuration 
is a sum of N single-particle states. The latter are eigen- 
states only for the free Hamiltonian without interactions. 
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At this point, numerical methods are usually employed 
to diagonalize the physical Hamiltonian which includes 
the Coulomb interaction as well as other effects ||] . 

Here, we will take a different approach. First, we 
note that any configuration \j\ , . . . , jpj] might contain 
so-called blocks of adjacent angular orbitals which are 
all occupied. Such a block jiiji + 1> • • • >3i + n ~ 1 IS 
uniquely characterized by two numbers, namely the num- 
ber n of adjacent occupied orbitals, and the number J = 
■ .+ji+n—l)/n = ji + ^(n—l). It is clear that 
we can enumerate all the configurations also by the blocks 
they decompose into, [(m, Ji), (n 2 , J2), (nt,je)]. Note 
that J might be half-integer, and that (1,J) denotes an 
orbital which does not have an occupied neighbouring or- 
bital. Of course, 1 < I < N with the extreme cases given 
by either one complete block or only isolated single or- 
bitals. Since the angular momentum orbital also defines a 
radial quantization, electrons in a block are also radially 
adjacent. In some abuse of language we might view each 
such configuration, projected on the radial component, 
as a configuration of a spin chain where spins can either 
be occupied or not. This leads us to the one hypothesis 
we will make in this paper: 

Hypothesis: The electrons within a block (n,j) are 
(sufficiently strongly) coupled due to spin-spin and spin- 
orbit interactions to form a composite state of charge n 
and "spin" j, while there is no (or negligible) spin-spin 
and spin- orbit interaction among different blocks. 

Clearly, this hypothesis is the equivalent of a nearest- 
neighbour interaction approximation. The clue is that 
we can treat a block as a collective state. In particu- 
lar, antisymmetrizing is performed seperately for each 
block, which will influence certain multiplicities. Also, 
the Coulomb interaction between two blocks (rii, j\) and 
(^2, 32) is calculated classicaly as given by the potential 
between two rings of charge n% and n 2 with radii sfji and 
t/j2 respectively. Here and in the following we adopt a 
system of units in which all proportionality constants are 
put to one, since we are only concerned with qualitative 
features. However, in most computations only ratios ap- 
pear, making them independent of our particular choice 
of units. Our choice is motivated by keeping everything 
as close as possible to pure integers emphasizing the com- 
binatorical nature of our approach. 

Our basis of composite states does already account 
for spin-spin effects, the remaining other important ef- 
fect being the Coulomb interaction. However, to further 
simplify all calculations, we make the following obser- 
vation: It has been shown numerically that the quali- 
tative behavior of classical dots is quite independent of 
the precise nature of the repulsive interaction V(ri, r,-) = 
eiej\ri — £j| -7 , with 7=1 yielding the Coulomb case, 
see e.g. jaj. The same should hold for quantum dots, 
as our results confirm a posteriori for 7 = 2. This is 
called universality, a property which quantum dots share 
with Laughlins wave functions for the fractional quantum 



Hall effect at filling factors v = l/(2p+l). However, the 
formulae for the potential energy between blocks and, in 
particular, the self-energy of a block, can only be given 
in a simple closed form in n and j for 7 even. Since the 
spirit of this paper is simplicity, we will therefore put 
7 = 2 with the silent understanding of universality. We 
have checked for a wide range of N, J that the qualitative 
features are indeed not affected by this. 

Thus, within our setting, the self-energy of a block 
(n, j) is given, in the equilibrium state, by the potential 
energy of a regular n-gon of radius y/j, which yields 



( n ,j) 



(n + l)n(n — 1) 
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Similarly, the potential energy between two such blocks 
is simply 



E, 



(™1 Jl),("2 ,32) 



nin 2 
\ji ~ h\ 



(5) 



Note that the squared radii in these formulas have been 
replaced by the "spins" of the composite states, bringing 
everything down to rational numbers. It will be imporant 
later that, if we were to include the precise proportion- 
ality factor of the relation j = m e ujr 2 , it would yield a 
common factor to all these partial energies. We are now 
ready to proceed in finding magic numbers. 



II. MAGIC NUMBERS 

All we have to do is calculate the energies for all block 
configurations according to the above formulae, where we 
neglect any pure quantum mechanical effects (such as the 
rest energy ^Nh) as well as the kinetic term (which is 
linear in J) yielding the so-called excitation energy 

E m ee E[(mJi), (ntJt)] 
t 

= ^2 E (m,3i) + ^2 E (n z J t ),(n kJk ) , (6) 
i— 1 i<k 

where we have assumed that we have enumerated all 
configurations (in a completely arbitrary manner) by m. 
We will now assume that the normalized ground state of 
the quantum dot can be written as a linear combination 

l*> = E™ S ™ K P(-l^m)l*m). with the I*™) forming 

an ortho-normal system, i.e. each configuration m will be 
weighted accordingly by its Boltzman factor exp(— f3E m ). 
The inverse temperature (3 is as yet unknown, but can 
easily be determined by solving the normalization condi- 
tion 



Z(/?)=^S m exp(-/?£ m ) = l 



(7) 



Please note that this is much easier than it looks! Putting 
x = cxp(— (3), one only has to find the unique real zero 
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of the polynomial Z(x') — 1 which lies in the interval 
x' G [0,1]. Here x' = x x l M with M the least common 
multiple of all denominators of the rational(!) numbers 
E m , and the polynomial Z(x') — 1 £ Z[x']. Finding this 
one zero can be done to arbitrary percision and very 
quickly, even for huge polynomials. The S m are sym- 
metry factors which take into account that each block 
is antisymmetrizcd individually, but also that each block 
(n, j) has a discrete Z„ symmetry. Hence individual an- 
tisymmetrization yields for each such block a symmetry 
factor n\\L n \~ 1 = (n— 1)!. This has to be multiplied by 
the number of ways to distribute the TV electrons among 
the £ blocks of a given configuration. Hence 



S m = S[(ni,ji), . . . , (nt,jt)] = 



N\ 1 

jn^- 1 ) 1 ' ^ 



such that the probability of the m-th configuration is 
Pm = 5 m exp(-/3S m ). 

Alternatively, we could leave (3 undetermined, treating 
it as an outer parameter of the system. This amounts in 
allowing different values for the normalization constant 
Z(j3). Again, we have checked that the qualitative fea- 
tures vary weakly with (3. Fixing (3 as above yields - in 
a loose sense - the quantum mechanical (inverse) tem- 
perature, a measure of how likely the system fluctuates 



quantum mechanically between its different states (ne- 
glecting thermal fluctuations due to a heat bath). 

After determining (3, we have the complete parti- 
tion function of our quantum dot system which we now 
treat as a system in statistical mechanics. Therefore 
we can easily determine observables such as the en- 
ergy (E) = dplog(Z) = J2 m PmE m , the specific heat 
(C) = I3- 2 C V = 9|log(Z) = (E 2 ) - (EY 2 , and the en- 
tropy S = - J2 m Pm log(p m ). As figures |l| and || demon- 
strate, the magic numbers show up in these observables 
exactly as they are supposed to do. We compare the most 
interesting cases, TV = 5 and N = 6. The temperature 
1/(3 and the specific heat are multiplied by suitable fac- 
tors to improve their appearance in the plots, ft should 
be noted that all these calculations, whose results are col- 
lected in the plots, only take a few minutes with a small 
and simple C-program on a typical workstation. 

Naturally, the "smoothness" of the curves shown in 
figures [l] and || increases with TV and for each TV with 
J due to the fast growth of cn,j- However, we find 
it noteworthy that a mesoscopic system of only 2 or 
3 particles already shows the same features, i.e. can 
successfully be approximated by a statistical system. 
One can interpret this fact that quantum mechanics in a 
mesoscopic system shows itself mainly by moving the sys- 



Thermodynamic data for quantum dot with N=5, J=10..60 



Thermodynamic data for quantum dot with N=6, J=15..75 
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FIG. 1. Thermodynamic data for a parabolic confining 
quantum dot with N — 5 electrons and total angular mo- 
mentum 10 < J < 60. The magic numbers at which the 
energy as well as the temperature exhibit local minima can 
clearly be seen. At these numbers, the entropy also has lo- 
cal minima indicating that the system has higher order for 
these angular momenta. Indeed, the system is closest to a 
cristalline structure at the magic numbers. The behaviour of 
the entropy is contrasted with the logarithm of the number 
of configurations. 
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FIG. 2. Thermodynamic data for a parabolic confining 
quantum dot with iV = 6 electrons and total angular mo- 
mentum f 5 < J < 75. Compare the magic numbers with the 
ones for TV = 5. A close look reveals that in addition to the 
period TV = 6 of magic numbers, equation (^), another period 
of TV — 1 = 5 is superimposed according to [3] which coin- 
cides with the former for J = 45. Of course, the same is true 
for the N = 5 quantum dot, but less prominently expressed. 
Note that magic numbers can already be seen for very small 
J ~ Jmin = TV(TV - l)/2. 
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tern from a purely classical configuration to a statisti- 
cal one. To show more clearly what we mean by this, 
we will consider the radial probability distribution, i.e. 
p(R)n,j = / o 27r |*(i?,0K,/| 2 d</> for given N, J. The 
point is that we do not even have to calculate anything 
of the wavefunction ^(i?, </>)n j. Since angular momen- 
tum and orbital radius are interlinked, we only need to 
consider the formal sum 

%) = E^fi^ i = X>^> ( 9 ) 

m i—1 j—Q 

where again m stands for the m-th configuration of 
blocks [(ni, Ji), . . . , (ri£, Ji m )\- The mathematically in- 
clined reader may wish to consider this as some kind of 
discrete Fourier transform of the radial or, more precisely, 
angular momentum probability distribution. The radial 
distribution can then be read off as p{R = y/j) = pj. 
However, we assumed here that the composite block state 
(n, j) is localised at j. Since this state is a collective state 
out of n electrons occupying the orbits from j — to 
j+ ^^-j it will certainly be smeared out over this region 
in angular momentum space. In the spirit of treating 
the system as a statistical mechanics ensemble, the most 
natural distribution of the composite state is a binomi- 
nal distribution centered at j with width ^ . One could 
equally take the point of view that a Gaussian distri- 
bution is more natural for a quantum mechanical state. 
However, since we are only interested in a qualitative 
analysis, we may neglect the difference between these two 
distributions (one may convince onself that both distri- 
butions lead to very similar functions p(r)). Since j may 
be half integer, the correct binominal distribution yields 

m i=l l i= v ' 
J 

= 5>y (io) 

as the radial distribution function. We are now ready to 
discuss the issue of symmetry loss, which is observed at 
the special magic numbers eq. (||) for N > 6. 

III. SYMMETRY LOSS 

The smallest non-trivial special magic number for N = 
5 is J(5,4) = 30, and for N = 6 it is J(6,5) = 45. Di- 
rect numerical analysis shows that around these numbers, 
(N — l)-gon symmetry with an occupied center is the 
preferred symmetry, whereas at the precise special magic 
numbers, and N > 6, symmetry is even more weakly 
expressed. One speaks of a fluid-like state at these num- 
bers (in analogy to Laughlins quantum droplets), and 



says that the quantum dot with these particular values 
of total angular momentum J = (N, k'(N— 1)) has filling 
factor 

= I N(N-l) _ ^N(N-l) _ 1 

2 J N(N-l){k'+±) 2k' + 1 ' 1 ' 

Before studying the radial distributions, we would like 
to explain this fact, and also that it is only observed for 
N > 5. For smaller particle numbers, Zjv symmetry is 
not broken (or only very slightly). Our considerations 
are made from very simple general assumptions and do 
not rely on any numerical studies. The key ingredient is 
the following relation between different magic numbers: 

J{N,(N-l)k')= (12) 
J(L, (L - l)k') + J{N -L,{N + L- l)k' + L) . 

Note that this relation can be applied recursively! The 
meaning is that for the special magic numbers, there are 
several ways to build a state out of smaller magic units. 
However, the above relation is quite different from the 
trivial relation J(N, k) = J(L, k) + J(N~L, k+L), which 
just divides one block of angular momentum eigenstates 
into two adjacent ones. The former relation is to be un- 
derstood as follows. If the offset k is divisible by N — 1 
and sufficient large, Z^v symmetry can be replaced by 
symmetry (with a rather small offset) within Zjv-l 
symmetry (with an even larger offset). This is possi- 
ble, because this distribution of angular momenta among 
the particles gives the same magic number, but with two 
sufficiently separated blocks to allow separate antisym- 
metrization. 

Clearly, this alone does not force loss of symmetry. 
The question is, which of these possible distributions has 
the lowest energy, and how close do the energies of the 
other distributions come to the minimal one. To an- 
swer these questions, one only needs the formulae (Q), 
, and (|^), and compare energies of different configura- 
tions. If the system were classical, we would only have to 
compare E[(N, J/N)} with E[(N - 1, J/(N ~ 1)), (1, 0)]. 
The former energy is smaller than the latter as long as 
N < i(13 + V73) < 6. Since a quantum dot is quantum 
mechanical by definition, all possible configurations con- 
tribute and we have to calculate the full partition func- 
tion Z instead. 

The above formula (|l2] ) is recursive. Moreover, one of 
the magic numbers on the right hand side is again a spe- 
cial magic number. Thus, it appears that the fluid- like 
quantum state is due to the fact that (for "large" N) 
there are a lot of different antisymmetrizations possible 
at these numbers which - as we are going to show - have 
all comparable similar energies, where the energy of the 
center occupied (N — l)-gon configuration is the minimal 
one for N > 5. There even exist antisymmetrizations 
with three or more nested blocks, whenever N is large 
enough. 
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We can read off the loss of symmetry from the radial 
distribution function. At a magic number J, we would 
expect it to have a well defined peak centered at R = \J~J 
for N < 5, and for N > 6 it should show two equally high 
peaks, one at R = 0. However, we find something quite 
different from that. In the following plot (figure ||), we 
compare the radial distributions for a range 39 < J < 51 
around the special magical value J = 45 for N = 6. 
One clearly sees that at the ordinary magic numbers 
J = 39, 51 full T,§ symmetry is restored (since there is 
no peak at R — 0)! Z5 symmetry with one particle in the 
center is most prominently exhibited for the neighbour- 
ing values J — 40, 50. This very sensitive dependency on 
the value of J is a pure quantum effect contradicting clas- 
sical expectations. However, the behaviour at the special 
magic number J — 45 is different. There is a certain 
probability for a particle to sit in the center, but it is 
only half of the probability of the main peak. This tells 
us that the main peak cannot come alone from the con- 
figuration with a five particle block and one particle in 
the center. 

One way to interpret this result is to say that a quan- 
tum dot behaves either more as a quantum mechanically 
system or more as a statistial ensemble depending on the 
value of the filling factor v defined in (|ll"| ) . The system is 
very much like a statistical ensemble and most resembles 



a quantum droplet for v — 2 w+i w ith £ Z + which pre- 
cisely happens for J being a special magic number. That 
is the same condition as for the filling factor for Laughlin- 
type fractional quantum Hall states which are described 
via the one-component plasma analogy as statistical sys- 
tems! The system is very quantum mechanical, and in 
fact resembles most an atom-like structure, when J is a 
non-special magic number. This does not translate into 
a straightforward statement with respect to the filling 
factor. 

The best way to compare the behavior at special magic 
numbers (||) with the one at ordinary magic numbers (|l]) 
would be to calculate the 2-point charge density correla- 
tion p(r, ro). In numerical studies of quantum dots, one 
usually freezes the position of one electron at ro and com- 
putes the probability of finding a charge at r. If ro is cho- 
sen such that it lies in an angular momentum orbital of 
high probability, one will get plots which clearly exhibit 
the symmetry properties of the quantum dot state. How- 
ever, to do this we would need the wave function which we 
avoided calculating. But the information encoded in ( |l0| ) 
is sufficient to give us a semi-classical approximation of 
the charge distribution. We simply pick the configuration 
with the lowest energy and put us into its "rest frame" . 
By this we mean the following: The minimal energy con- 
figuration has a certain block decomposition, and for this 



Smoothed radial probability distribution p(R) for N=6, J=39..51 



or 




FIG. 3. Some radial distributions for a quantum dot with 
TV = 6 electrons and total angular momentum 39 < J < 51. 
Magic numbers are J = 39, 45, 51 with J = 45 being a special 
magic number. Note that the classically expected 5-fold sym- 
metry is only observed for J = J mag ± 1. This is because the 
system prefers forming composite states which works best for 
magic numbers. However, the neighbouring J-values behave 
similar to atoms with one shell either almost filled or filled 
plus one excess electron. 




FIG. 4. Comparsion of approximated charge distributions 
(in the rest frame of the minimal energy configuration) for 
N = 5 (top) and N = 6 (bottom). Left are special magic 
numbers corresponding to v = 1/3, right are ordinary magic 
numbers. Note the difference for N = 5 and N = 6, the latter 
showing much stronger quantum droplet behavior at a special 
magic number. 
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decomposition we localize the electrons on the edges of 
the appropriate n-gons, smeared out according to the bi- 
nominal distribution taking into account the quantum 
mechanical width of the n-gons. For TV < 5 the minimal 
configuration is simply one AT-gon, for 6 < N < 8 it is 
a (N — l)-gon with an occupied center. For N > 9 we 
have a shell like structure where the relative orientation 
of the electrons in one shell to the ones in another would 
have to be chosen according to the classical solution. All 
other configurations are considered to be radially sym- 
metric. In this way, superimposing these approximations 
of charge distributions, we obtain a measure of how much 
the minimal energy configuration sticks out of the other 
ones, and hence a measure of how much the particular 
symmetry of the minimal energy configurations deter- 
mines the symmetry of the full charge distribution. When 
plotting charge distributions in this manner, it is clear 
that the integrated radial probability for a block within 
a configuration (first equation in |Io| ) has to be weighted 
not only by the binominal distribution, but also by l/2irr 
where the radius r varies over the width where the block 
is localized. Moreover, such a charge distribution is natu- 
rally time averaged. To be more precise, we approximate 
the probability distributions as follows: A block (n,j) 
with probability p, if it does not have the dominant sym- 
metry of the minimal energy configuration, contributes 



1 



(13) 



where T> a (x) denotes a symmetric distribution centered 
around zero with width a, in our case a normalized 
binominal distribution such as V a (x) = 2~ 2o T(2(7 + 
l)/r(tr + I — x)T(a + 1 + x). Since this is independent 
of cj), it is radially symmetric. Of course, as mentioned 
above, we could equally well have chosen another distri- 
bution such as a Gaussian one without affecting the final 
results very much. In the case where the block has the 
dominant symmetry of the minimal energy configuration, 
we use instead 



PU> <l>) = P Vn / 2 (Vj 2 +f - 2jjcos(27r(fc/n- <f>)) 



k=0 



2VI.7-J1 + 1 



(14) 



Our approximation differs from the approach where one 
electron is frozen on a likely position only in the fact that 
we freeze the position of one suitable composite particle 
state. Therefore, our scheme is essentially only the conse- 
quent translation of computing pair correlations p(r, r ) 
in a single particle basis to our basis of composite particle 
states. As figure [| clearly demonstrates, full Af-gon sym- 
metry holds even for A^ > 6 for ordinary magic numbers, 
while symmetry loss occurs for special magic numbers 
and N > 6. 



IV. CONCLUSION 



Although space permits only to present these few ex- 
amples, we have performed calculations of quantum dot 
states for 2 < A^ < 9 and with total angular momentum 
up to J = 126. We can confirm the occurence of special 
magic numbers for N = 6, 7, 8 corresponding to filling 
factors v = 1/3,1/5,1/7, surpassing the range of N, J 
accessible to exact diagonalization methods. We are also 
able to see that for A" = 9 more complicated shell-like 
structures appear, but that magic numbers still lead to 
an enhancement of symmetry, and special magic num- 
bers to a weakening of symmetry. The precise results 
will appear elsewhere. 

We proposed a new and extremely simple way to un- 
derstand the qualitative features of quantum dots by mix- 
ing classical and quantum mechanical points of view. We 
think that the results justify our approximation scheme, 
where we chose a basis of composite particle states for 
which only the classical repulsive interaction is taken into 
account. Quantum mechanics entered the game mainly 
through the quantization of angular momentum which 
made the computation of the partition function a finite 
(and in fact quite simple and fast) enterprise. As a fur- 
ther simplification, we made use of universality and re- 
placed the Coulomb interaction by a V = \r\~ 2 potential. 
This allowed us to perform all computations with integer 
arithmetics alone. Once one accepts our special choice of 
a basis of states, all further computations can hence be 
done exactly within the scheme. Resorting to the physi- 
cal Coulomb potential makes computations slighly more 
complicated but does not, as we have checked, change 
the qualitative features of the results. The main influ- 
ence of varying 7 in V = |r| 1 is that the relative radial 
positions of blocks are shifted somewhat. More precise 
statements about universality can be found in ||. 

Of course, our approach is over-simplifying and cannot 
capture finer details of quantum dots. There is one par- 
ticular point of possible criticism: We assumed that par- 
icles within a block would distribute themselves equidis- 
tantly on a circle. Classically this is not true in general if 
the configuration consists of more than one non-trivial 
block. However, our central hypothesis suggests that 
quantum mechanically each block should be treated as a 
smooth delocalized ring shaped charge distribution, since 
we assume that each block forms a collective state. In- 
deed, when calculating the potential energy between two 
blocks via (|J), we assume precisely that. On the other 
hand, when we calculate the self-energy of a block, we re- 
sort to a classical picture of n point-sources distributed 
equidistant on a circle, see (0). It is here, where we are 
in danger of a systematic error, slightly underestimating 
the self-energy in the presence of another block (which 
may force the electrons to move to a non-equidistant dis- 
tribution). However, if our central hypothesis is true, it 
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depends on the strength of the coupling of the electrons 
within a block how much the resulting collective state 
is influenced by the presence of another block. In our 
simplistic picture, we assume that such corrections are 
neglegible and do not change the qualitative patterns. 
This is plausible since spin-spin interactions among com- 
pletely polarized electrons lead to a further repulsive in- 
teraction which locally, within a block, might dominate 
effects coming from the electro-static interaction with 
other blocks. 

Another possible point of criticism is our normalization 
condition ([?]) which assumes ortho-normality of our basis 
of states. This is of course not entirely true. However, 
figures 1 and 2 clearly show that the resulting temper- 
ature estimates are not too far off what we would ex- 
pect, namely (E) ~ NkT for iV particles. In fact, one 
can easily replace (0) by the condition Z{[3) — C/3~ N 
with C chosen such that for the minimal possible case, 
J = N(N - l)/2, Z(J3) = 1. Solving this constraint 
is slightly more complicated, but the other results differ 
not very much from the ones obtained with our idealistic 
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FIG. 5. Comparsion of approximated charge distributions 
(in the rest frame of the minimal energy configuration) for 
N = 6. Left are the distributions for filling factors v — 1/5 
(top) and v = 1/7 (bottom) with J = 75 and J = 105 respec- 
tively. This is compared with the distributions for non-special 
magic numbers close to the special magic ones on the right. 
Note that no full restoration of symmetry occurs for J — 81 
(top), where the j — 1 orbital can be seen to be occupied, 
but where the lowest-energy configuration barely sticks out 
of the more likely fully symmetric configuration. In case of 
J = 99 (bottom), full symmetry restoration is still observed. 
See main text for details about which magic numbers fail to 
completely restore maximal symmetry. 



approach. In particular, the radial distributions just be- 
come a bit more smeared out, but retain their qualitative 
structure. 

Our results bring us to the conclusion that our hypoth- 
esis is a good one, i.e. that it predicts the correct qual- 
itative behavior expected from quantum dots. It would 
hence be very interesting to derive this hypothesis from 
a microscopical treatment of quantum dots. It is well 
known that spin-spin and spin-orbit interactions are im- 
portant in the theoretical understanding of electron or- 
bitals of atoms. Quantum dots are often described as 
"artificial atoms" ||, and it would only further justify 
this labeling, if they indeed shared the importance of 
spin-interaction effects with their name-inducing natural 
relatives. Also, it would be interesting to relax our con- 
dition of completely polarized electrons in order to study 
quantum dots where angular momentum orbits can be 
occupied by upto two electrons. Finally, a more realis- 
tic normalization condition for the partition function, as 
mentioned above, might be desirable. But these investi- 
gations will be left for future work. 
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FIG. 6. Comparsion of approximated charge distributions 
for different rest frames. Top row shows the quantum dot for 
TV = 6 and J = 81 in the rest frames for full 6-fold symmetry 
(left) and for 4- fold symmetry (right), to be compared with 
the plot for 5-fold symmetry in figure 5 (top right). The fully 
symmetric configuration is still the dominant one. Bottom 
row shows a similar comparsion for N = 7 and J — 70 (the 
first magic number after the special magic number J = 63 
for v = 1/3). The difference between the lowest energy rest 
frame (left) and the highest symmetry rest frame (right) is 
very small, but the latter is the dominant one. Hence, symme- 
try restoration still takes place, but is much weaker expressed 
than for other non-special magic numbers. 
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V. PREDICTIONS 

So far, we have successfully reproduced results which 
were already achieved by other methods, mainly by exact 
diagonalization techniques. Since our method is simple 
and fast, we can probe higher particle numbers and to- 
tal angular momenta than the ones accessible to numer- 
ical methods. We will present here some calculations for 
-/V = 6 with much larger total angular momentum as can 
be found in the literature (see figure ||) . 

In particular, we confirm that symmetry loss - as ex- 
pected - occurs for higher special magic numbers cor- 
responding to fillings v — 1/5 and v = 1/7. On the 
other hand, full symmetry is still restored for most but 
not all non-special magic numbers! This is a new phe- 
nomenon, and also a nice demonstration of mesoscopic 
physics. Most non-special magic numbers do not con- 
tain states of the form [(1, Ji), (N — FJ2)] with energies 
comparable to the (for magic numbers always existing) 
configuration [(N,J — J/N)). The reason is that for 
most magic numbers there is no choice with very small 
ji which would favor a smaller energy than the energy of 
the [(N,J — J/N)] configuration. One of the properties 
which make special magic numbers special is that at these 
values of the total angular momentum the configuration 
[(1, 0), (N - 1, J/(N - 1))] becomes possible. 

However, if the total angular momentum J is suffi- 
ciently large, configurations such as [(1, 1), (N — 1, (J — 
l)/(iV — 1))] may have energies smaller than the en- 
ergy of the most symmetric solution. In the case of 
N = 6 this happens the first time for J = 81, where 
the above given configuration consisting of a block with 
5 particles around a single particle orbiting the center 
in the j = 1 orbital has indeed a slightly lower en- 
ergy than the configuration with a full block out of 6 
particles, although the former state has a much smaller 
probability than the latter. This is simply due to the 
fact that for J ^ 1, the difference between the ener- 
gies of the configurations [(1, 0), (N — 1, J/(N — 1))] and 
[(1, 1), (iV-1, (J-l)/(iV-l))] becomes arbitrarily small. 
Put differently, for J ^> 1 the system behaves more 
and more classical, but the transition from the quan- 
tum realm to the classical world is not continuous with 
increasing J. For example, full "Lq symmetry is restored 
for the magic number J = 99. More precisely, the first 
non-special magic numbers (Q) where symmetry restora- 
tion does not take place are J(N,k'N) for k' > k' Q (N). 
We have for example k' Q (6) — 2 and k' (7) = 1. These 
happen to be the non-special magic numbers which fol- 
low next to special magic numbers. For even higher J we 
might expect that more and more magic numbers fail to 
restore the full symmetry of the system. 

Because the minimal energy configuration does not 
possess the maximal possible symmetry, it may have a 
smaller probability than the maximal symmetric config- 



uration due to different multiplicities. In fact, this is 
the case when the magic numbers are non-special. When 
the magic number is special, the minimal energy con- 
figuration is the dominant one, although the dominance 
is not strongly expressed - which is precisely the reason 
which leads to the droplet like probability distributions 
of such quantum dot states. It follows, that symmetry 
restoration should still work to some degree for the above 
discussed non-special magic numbers. We demonstrate 
this in figure |6|, where we compare probability distribu- 
tions for different rest frames, namely the rest frame of 
the minimal energy configuration and the one for the 
maximal symmetry configuration. One can infer from 
these plots that symmetry restoration does take place 
for the non-special magic numbers J(N,k' N), but of a 
somewhat lesser degree than for other non-special magic 
numbers. 

Therefore, we have demonstrated that our method is 
capable of probing unknown regions in "quantum dot 
space" showing and predicting new patterns. We be- 
lieve that our method is a useful tool to explore in more 
detail how quantum mechanics and classical physics are 
interlinked in mesoscopic systems. 
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